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Abstract. 


Seismic records from IRIS (NIL, KBL, MSEY, KMI, QIZ, DGAR) were used to 
compute AWARAN September 24 2013 Earthquake's focal mechanism. Earthquakes are 
one of the most drastic threats that are likely to cause heavy human casualties and can 
destroy entire cities within minutes. We processed the waveform data from 6 separate 
stations for the computation of displacement. The issues related to attenuation have been 
overcome by the use of low frequencies (0.05-0.1). The crustacean model that determinesthe 
tensor PREM model's time is determined by minimizing the variance but produces a robust 
focal mechanism. Inverting the displacement data is applied to recover the moment tensor 
from the source, calculated the green feature with the reversal. GIS-based vectors offault lines 
were abstracted by the digitization of regional tectonic lineaments and structuralmaps. We 
developed synthetic and observed surface shapes with a 100 percent dc correlation. Our 
understanding of the features of crustacean deformation and regional tectonics in the area 
will benefit from the resulting focus mechanisms. The measured earthquake focal 
mechanism is similar to a USGS agency hit, dip, and rake. 

Keywords: Awaran earthquake, Focal Mechanism, PREM Model, Green Function, GIS. 
INTRODUCTION 


Seismic tremor or earthquakes are natural hazards which devastate the natural beauty and 
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ecological balance of a region badly, [1]. The earthquake in AWARAN on September 24, 2013, 
was ttiggered by an oblique landslide movement in one of Chaman 's southern strands, [2]. The 
event took place within the large area of the plate limit established by the MakranCumulative 
Prism and the Chaman Fault System, [3]. It was situated in a region with well- 
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defined fault morphologies, [4 ][5 ]. but large (>M7) historical earthquakes were not 
previously observed, | 6 ][ 7 ].The surface crust depth (15 to 20 km), the northern slip 
movement and the internal depth was approx. 80 km as shown [8]. In 1935, the Mw 7.5 
earthquake occurred at the Ghzaband fault, [9]. The core mechanisms of medium-depth 
earthquakes have been clarified mainly in response to the decreasing extension in many 
subduction areas worldwide, [10]. Seismic procedures provide the highest accuracy of the 
earth's internal structure than any geophysical method, [11]. This is because elastic waves of any 
geophysical wave have the shortest wavelength and because the physics which regulates their 
spatial and temporal sensitivity determines the origin of energy, [12][13]. 


Several methods for evaluation of point-source and source-time (STF) process for 
events were established simultaneously, [14]. The use of waveform modelling for calculating 
focus mechanisms depends on modern tools [15][16][17]. The tensors depend on the spring 
strength and fault path [18]. The pivoting function or the solving of problems of the plane 
also implies "beach balls." It gives us a feeling about seismic waves that originate from a 
complex rupture, involving space and time displacements that are abnormal [19]. 
Furthermore, it is a plus point to write that is any orthogonal coordinate system since vector 
and tensor equations are valid irrespective of the coordinates of system. Ordinary data 
reversal techniques [20] surface waves [21][22][23] and body waves [24][25] have been 
evaluated and the reversal of broadband seismic charts has been proposed for regional 
earthquakes [26]. Other methods are also used to evaluate influential movement records and the 


depth of the source [27]. 


The ISOLA Interface is used for instance, which carry out certain tasks to access all 
investment procedures, [28] [29]. The green function is measured by wave frequency 
method, [31] [32]. The code is written in Fortran, but the user can control its use through 
ISOLA-GUI, [33]. ISOLA-GUI creates all the required files with the inspection diagrams 
(Descartes and Geographies coordinates) for calculating the Green-Inversion function, [34] 
[35]. The relationship representing the percentage of a double pair solution is defined, since it 
can show problem output in investment. The modules run concurrently and most of them 
"know" the answer from the user or try to suggest acceptable processing parameters. 
Significant information on regional and local tectonics is given for medium sized earthquakes in 
the country, [36]. The focusing mechanisms calculated using only the first direction of 
movement can thus have a significant effect on the parameters of the coordinating 
mechanism calculated by these incorrect measurements of the first movement. Depending 
upon the distribution and consistency of the first motion data, more than one focus 
mechanism fits to the data. The key benefit of the instantaneous tensor is that seismographs 
can be reversed to identify the source parameters. Instant tensors for reflection are definedas 
the movement of an error during an earthquake consisting of nine pairs or nine sets of two 
vectors. 

Material and Methods. 
Investigation site. 


Sep 2020 | Vol 2|Issue 3 Page | 110 


oven (access International Journal of Innovations in Science & Technology 


There were many earthquakes, mostly occurred in coastal areas in Balochistan. District of 
Awaran is situated (Fig. 1) in the south of the province of Baluchistan (Lat: 27,016 ° N; 
65,547 ° E). The red zone for earthquakes in Pakistan is known to be the Mekran range dueto 
tectonic activities. In the remote and inaccessible regions, the earthquake caused widespread 
destruction and the death toll reached over 825 . Awaran, Ritaj, Jashkur, Nok Jo, Parwar, Dandar 
and Hoshab were the affected areas. Almost 33,000 homes and 300,000 individuals suffered 
from the earthquake, [37]. . 

The Awaran earthquake was caused by the southern strands of the Hoshab fault in 
the Chaman fault system. An analysis of high-resolution World View 2 satellite images also 
confirms the rupture on the left side of the surface. Another confirmation was made of 
Landsat-8 image for the rupture that corresponded to the Hoshab failure. 

The event ruptured the Makran subduction megathrust. The earthquake was 7.7 as 
measured on the Richter scale. A new island formed called “mud volcano” on the sea floor 
following the earthquake. The epicenter was situated in an area influenced by the above two 
crash logics. A complex collection of mountains folding and fissured along the coast and 
northbound to the west from India and Eurasia's convergence area resulted from relative 
closeness between the two species [38]. 
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Figure 1. Location of Awaran earthquake September 24", 2013Awaran 
Earthquake in Terms of Richter Scale 


Charles F. Richter developed the Richter Scale in 1935. Originally, the earthquake was 7.4MW 
and then upgraded to 7.8 but later changed to 7.7, on the Richter scale. A heavy shaking 
(VII) occurred in the nearby town of Awaran on the 24th September, with a major shakingin 
the big town of Kharan and a heavy shaking in Hyderabad and mild shaking in Karachi. 


Table 1. Details of September 24 earthquake 
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Moment Magnitude Depth Percent Half Catalog Data Contributor 
of DC Duration Source 


4.418e+20 7.7 Mww 23.0 96 % US US? US? 
N-m km - 


Nodal Planes 
Table 2. Nodal Planes of the Study Area 


Plane Strike Dip Rake 
NP1 228° 46° 6° 
NP2 134° 86° 136° 
Principal Axes 
Table 3. The Details of Principal Axes of the Area 
Axis Value Plunge Azimuth 
T 4.458e+20 N-m 33° 81° 
N -0.080e+20 N-m 45° a9" 
P -4.377e+20 N-m 26° 189° 


Tectonics of the Area 


Makran tectonic map is shown in the figure 2. 
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Figure 2. There are major faults on the tectonic map of southern Pakistan. A full 
black circle is the core location of the Awaran earthquake [18] 

It spreads over 400 km north and a 1000 km east, and can subducts under the 
Eurosian platform by 3 cm / year north of the Arab plate and Ormara.. In 1945, an 8.0 MW 
eatthquake broke the enormous Makran thrust. The Chaman fault system is a 1200 km left- 
hand failure system adjusted for 3 cm / yr north of India in relation to Eurasia. The Chaman 
fault system consists of the Ornach Nal, Ghazaband and Chaman 's sinister failures. In the 
south, the Chaman defects refer to the Seyhan, Punjur and Hushab Makran chains! 


Sep 2020 | Vol 2|Issue 3 Page | 112 


oven (access International Journal of Innovations in Science & Technology 


defections and those defections are internalised; the tale from India and Eurasia and the 
shorter one from the Arab Peninsula and Eurasia. The faults left in the side blow sliding 
faults of the Chaman and Urnash-Nal power the Indians and Eurasians' crash. Landslide 
systems predominantly form the Makran subduction zone boundaries. The fault system 
Ornach-Nal marks the eastern frontier, while the fault system Minab-Zendan represents the 
western border. As the Arabian plate is compressed, in the Makran area and the Chagai Arc, 
there is a directional tendency between east and west. It implants with an implantation angleof 
2-8 degrees at a speed of 19.5 mm pet year. 


Figure 3. 2013 Earthquake in MW 7.7 Awaran. (a) An aerial view of the Oran area's 
destruction. (b) A view of the island following an earthquake in Awaran [8] 

In Ziarat, the 1.300 Kilometers long Hass Soliman fracture, the Chaman fault and 
the Sibi basin are regulated by the tectonic and structural patterns. Quetta's local defect 
structures are comprised of Chaman 's Fault, Chengin Gulat structure, Harnai Tarta 
Structure, Gzband Gob Chassis, and Mach Chassis. The Makran Coast Fault is 225 km longin 
south of Ormara and Pasni. The earthquake occurred in the depth of the ground crust over 
the Makran subduction zone as a result of oblique displacement. 

The study is based on the seismic waveforms data which is available at 
https://ds.iris.edu/ds/nodes/dmc/data/types/waveform-data/ that was recorded at 
different stations of the world in response to the earthquake that struck just near the district 
Awaran on September 24, 2013. 


Table 4. Stations Information of September 24, 2013, AWARAN earthquake 


Station name Latitude Longitude 
NIL 33.65 73.269 
KBL 34.541 69.043 

DGAR -7.412 72,452 
MSEY -4.674 55.479 
QIZ 19.029 109.844 
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KMI 95 193 102.74 


For proper coverage of Awaran earthquake, we set the distance from 0 to 180 
degrees using powered. The coverage to the event was ensured from all sides for a better 
display of results. The start and end time are adjusted such that complete waveforms are 
obtained. A total 6 IRIS stations for Awaran earthquake were selected such that a proper 
coverage to the earthquake epicenter was ensured. Our goal was to measure the tensor 
moment. The first-order seismic tensor is an extensive representation of earthquake sources 
resulting in a straightforward linear relation with a sequence of green and fundamental error 
reactions. . So, we choose the station for which waveform data was available for all N, E, Z 
and time components. 
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- Figure 4. In the figure the orange circle shows the location of September 24th,2013 
Awaran earthquake epicentre, While the red triangle shows the location of the 
tecording NIL station (located in Nilore Islamabad) 
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Figure 5. Flow chart of Methodology Adopted for the work 


Defining Proper Crustal Model for Earth 


ISOLA-GUI operation is a crustal model built with the crustmod tool for all 
subsequent analysis steps (for example, green function measurement, polarity control, etc.). The 
vector speeds (Vp and V), density, P and S (Qp and Q) are inserted either in text boxesor in 
text files when the Matlab code provides correct formatting for ISOLA fortran coding.Some 
devices were available with experimental relationships for generating Vp and V velocity 
graphs with depth and measurement of V and strength. For each season, the symbol supports one 
cortical model, but future releases support several cortical models. We used the PREM Crustal 
model to calculate the current tensor for the study period. The main reference earth model 
(PREM) was 1-D used for seismic experiments on the planet for several yearsThis model has 
been developed for various data sets, including free centre oscillation frequency 
measurements, monitoring of surface waves scattering, travel times for different phases of body 
waves, tadio astronomical data from the field, mass and moment of tranquilly. For our study 
area, we used 12 layers of the earth's surface (depth 0) to the innercore (6371 km). 
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Figure 6. In this PREM Model above the red line shows Vs (Shear wave velocity) 
andblue line shows the Vp (Primary wave velocity). We notice that from 3000-5200 
the Shear wave velocity (Vs) value obtained is zero. 


Green Functions Computation 


The green function is computed using the waven number method. Like coupon 
delivery, several codes are available remotely to represent data in bullet source or multiple 
instantaneous tension models. Fortran codes exist as separate executable files called ISOLA- GUI 
for tasks such as computing and Green reversal. For additional M5 + earthquakes at 
distances up to 400 km, the code is also proven useful. The overall frequency for our workto 
measure the green function was 0.1. It was crucial for our results to use such a low frequency 
in our work, because the attenuation of such lower frequency is not very important. 
Therefore, the attenuation issue largely doesn't affect our performance. 

Spatial Analysis 


The primary data of earthquakes was accessed from open source earthquakedatabases 
of US Geological Survey and Pakistan Seismic Network of Pakistan Metrological 
Department.The dataset of 4031 events (spreading over 7 °N to 32° N and 55° E to 71° E) 
has been self-reviewed and annotated for cataloging. The magnitudes of those events were 
measured in different scales ranging between 3.7 <m< 8.1. All kinds of magnitude scales 
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have been translated into unified magnitude scale Le. moment magnitude (Mw) which is 
considered most reliable magnitude scale for cataloging. The capabilities of ArcMap tools 
have been explored for data analysis, clipping and mapping of desired earthquake data. The 
clipping tools of ArcMap were utilized to select the events within the buffer around a fault 
line and tectonic lineaments. The clipping tool improve the efficiency in getting data pointsin 
response of spatial queries. Since, the primary key of earthquake dataset is the spatial 


coordinates of epicentre of each event. 
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Figure 7: History of Earthquakes from 1900 to 2016 
Result and discussion. 


This research involves calculating the location and frequency of reflection of the 
waveform for a test depth in the horizontal source (in the SED position, i.e. in the USGS 
Focus centre). Fortran ISOLA code is used to revert the Waveform to use speed, while the 
user interacts even during reflecting in the Mat lab environment. The broad slip of a young 
failure and an active failure is a riddle about fault geometry's supposed impact on failure 
dynamics. The seismically active error plans tend to be optimally directed to Anderson's 
error theory. Slip or slip sources often come out from the axial mechanism of low-cross 
earthquakes, but slip slipping is relatively rare.. In Fig.8, which shows the stability of the 
focal length and the best acceptable depth at 17.5 km, the connexion between observations 
and artificial seismograms is shown in terms of depth. The maximum correlation was foundat 
0,43 at a depth of 17,5 km and the CC is 100% double. Baluchistan’s earthquake highlighted 
the engineering flaws, preconfigured stress and dynamic vulnerability, as was the case on several 
previous occasions on the seismic breakage propagation from dynamic modelling. 
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Correlation vs Depth Plot 
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Figure 8: As a function of the test source deepness down the SED focus centte, 
the correlation between the observed, artificial and focal mechanism waves. DC 
percentof colors 
We determined the earthquake focal mechanism, which equates the USGS strike, sinking and tilt. 
Then we reverse the seismogram in order to measure the time of launch, the intensity ofthe 
earthquake and the tensor's six components. Figures 9 and 10 demonstrate the effects ofthe 
reflection. Seismic moment was observed approximately 4,037e + 20 Nm, the equivalent of Mw 
7,7. US Geological Survey has reported the same scale. In contrast to a depth of 15.5 km, the 
best approach was measured at a depth of 17.5 km. The synthetic seismograms, asshown in 


the figure, well suit the observations. 
MOMENT TENSOR SOLUTION 
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Figure 9. Moment tensor solution of the Awaran earthquake. Triangles show the 
station used in the inversion. 
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Figure 10. Observed (green line) and synthetic (blue line) velocity (m/s): 
Frequency ranges 0.005—0.1 Hz, Event 1. Peak amplitudes (in m/s) are on the 
tight-hand side ofthe figure Since the occurrence of event, the times (in second) 
are given on the bottom of the plot 
Cross-section of Slip Distribution 


The strike direction is indicated above each fault plane and the hypocenter location is 
denoted by a star. Slip amplitude is shown in color and the motion direction of the hanging 
wall relative to the footwall (rake angle) is indicated with arrows. Contours show the rupture 
initiation time in seconds. 
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Figure 11. Cross-section of Slip Distribution 
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Surface projection 


Surface projection of the slip distribution was superimposed on GEBCO bathymetry. 
Thick white lines indicate major plate boundaries [Bird, 2003]. Gray circles, if present, are 
aftershock locations, sized by magnitude. 
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Figure 12. Surface Projection 
Spatial analysis of Epicentres of magnitudes 


The study area is seismically active, it is diversified in seismogenic potential, 
earthquake born processes, and consequential variation in magnitude strengths of the 
earthquakes. The magnitude of earthquakes quantifies the amount of elastic strain energy 
released during the earthquake. It was revealed that the changes in magnitude strength arenot 
just the statistical data variations, but these changes are closely articulated with the seismo- 
geological characteristics, deformation and depth of the active fault segments, microseismic 
proximity to the active margins, structural stability, fault mechanics and its kinematics etc. 


Sep 2020 | Vol 2|Issue 3 Page | 120 


oven (access International Journal of Innovations in Science & Technology 


GO'SOU"E E3400" GE*BD'O"E = GS*20'"E —72100"E —7SO'U"E = 77*S0'"E §—-BU*40'0"E 


soa" ersooL 


(lurkmenistan} 


28 0O'N 


31°50'0"N 34°40'0"N 37°30'0"N 
31°S50'0"N 34°40'0"N —-37°30'0"N 


29°0'0"N 


z 
2 
= 
R 


z 
z 
& 


"N  2610'0"N 
"N 26°14 0'O"N 


09°400°E en'soe 

Nottett| Weak | Ught [Moderate] Strong [Very strong 
pone | very oh] Light | Moderate 

a 9a | «2 | 2 | 2 rr a SE) 
7 


GOSOU"E ES*400"E GSSDDE 69'200"E —72100"E 7S'OO"E  77SDO'E  B0°400"E 


Figure 13. (a)Seismic map of Pakistan Figure 13.(b) Shakemap of Pakistan 


Tectonic Summary 


The Awaran earthquake of magnitude, M 7.7 earthquake in south-central Pakistan 
occurred as a result of oblique strike-slip type at shallow crustal depths. The location, focal 
mechanism solutions, and finite-fault modeling of the earthquake are consistent with left- 
lateral (southwest-striking) rupture within the Eurasia plate above the Makran subduction 
zone. The event occurred within the transition zone between northward subduction of the 
Arabia plate beneath the Eurasia plate and northward collision of the India plate with the 
Eurasia plate. The epicenter of the event is 69 km north of Awaran, Pakistan, and in 270 km 
north of Karachi. On a broad scale, the tectonics of southern and central Pakistan reflect a 
complex plate boundary where the Indian plate slides northward relative to the Eurasian 
plate in the east, and the Arabia plate subducts northward beneath the Eurasia plate in the 
Makran (western Pakistan). These motions typically result in north-south to northeast- 
southwest strike-slip motion at the location of the September 24th earthquake that is 
primarily accommodated on the Chaman fault, with the earthquake potentially occurring on 
one of the southernmost strands of this fault system. 
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Figure 14. Tectonic map of Pakistan 


While commonly plotted as points on maps, earthquakes of this size are more 
appropriately described as slip over a larger fault area. Strike-slip events are typically about 
160x20 km (length x width). Modelling of this earthquake implies dimensions of about 
200x20 km, predominantly up-dip and southwest of the hypocentre. Although seismically 
active, this portion of the Eurasia plate boundary region has not experienced large damaging 
earthquakes in the recent history. 


Conclusion 


The data from 6 stations (NIL, KBL, MSEY, KMI, QIZ, DGAR) was used to 
compute focal mechanisms. By using the low frequencies (0.05-0.1), the problems arisingdue 
to attenuation were significantly reduced. The maximum correlation was noticed at 0,43 at a 
depth of 17,5 km where the CC was 100% double. Baluchistan’s earthquake highlighted the 
engineering flaws, preconfigured stress and dynamic vulnerability, as it was the case onseveral 
previous occasions on the seismic breakage propagation from dynamic modelling. The plot of 
observed waveforms versus synthetic waveforms is well correlated. The use of a standard ground 
model such as PREM provided a robust focal mechanism. The Mat lab- based graphical user 
interface makes data management easy, while offering full user controlin the workflow and an 
intuitive graphical overview of the findings. The quantity of DC (double pair) was estimated 
as 7,7 100%. Today Pakistan is a seismically active and vulnerable to major seismic events. In 
Pakistan, tectonics needs to be re-invested as a matterof urgency and secret defects found in 
areas where the population is seriously endangered. 
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